Two-Dimensional Polymers with Random Short-Range Interactions 
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We use complete enumeration and Monte Carlo techniques to study two-dimensional self-avoiding 
polymer chains with quenched "charges" ±1. The interaction of charges at neighboring lattice sites 
is described by qiqj. We find that a polymer undergoes a collapse transition at a temperature Te, 
which decreases with increasing imbalance between charges. At the transition point, the dependence 
of the radius of gyration of the polymer on the number of monomers is characterized by an exponent 
vb = 0.60 ± 0.02, which is slightly larger than the similar exponent for homopolymers. We find no 
evidence of freezing at low temperatures. 
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A polymer in a solvent is subject to monomer- 
monomer interactions which consist of a short-range re- 
pulsion and a slightly longer range attraction. At high 
temperatures (or good solvent conditions), the repulsive 
interactions are dominant, and the polymer is swollen: 
Its radius of gyration Rg scales with the number of 
monomers N as Rg ^ N'^, with i/ > 1/2. As tempera- 
ture T is lowered (or solvent conditions worsen), a point 
is reached, called the 0-point, where the repulsive and 
attractive interactions effectively cancel and the polymer 
scales, for space dimensions d > 3, like an ideal ran- 
dom walk {v = 0.5) [0,1). For T smaller than the tran- 
sition temperature Tg, the attractive interactions prevail 
and the polymer collapses into a compact shape, with 
V = 1/d. Such homopolymers are often modeled by self- 
avoiding walks (SAWs) on a discrete lattice, with the at- 
tractive interactions included by introducing a negative 
energy for each pair of monomers residing on neighbor- 
ing lattice sites. Numerous Monte Carlo (MC) and exact 
enumeration studies of the 0-point have been performed 



The collapse transition of heteropolymers is particu- 
larly interesting in view of its possible relation to the 
problem of protein folding While models based 

on random heteropolymers significantly over-simplify the 
complexity of real proteins, they bring in fresh perspec- 
tives from the statistical mechanics of random systems 
and spin glasses [||J|] . A question of high interest (with 
no definite answer) is whether the interactions between 
different monomer types can modify the collapse transi- 
tion of a polymer 1 10|. Another interesting feature of het- 
eropolymer chains concerns the compact state: Is there a 
freezing transition at some temperature Ty(< Tg), below 
which the configurational entropy per monomer vanishes, 
and few conformations with low energy dominate [p|,pT|. 
This freezing transition is assumed to be analogous to 
the glass transition in Derrida's random energy model 
(REM) jlj], although this analogy is not always valid 

We study self-avoiding polymer chains on a two- 
dimensional (2D) square lattice. Each polymer chain 
is formed from two types of monomers, with charges 
Qi = ±1. The interaction Hamiltonian is: 



i,3 



(1) 



where = 1 if monomers i and j are located on ad- 
jacent lattice sites and = otherwise. The homo- 
geneous repulsion is imposed by the constraint of self- 
avoidance. The use of a 2D system enables us to inves- 
tigate longer chains than in 3D models and thus 
examine aspects not studied before. We find that a poly- 
mer undergoes a 0-transition, with a critical exponent 
i^e — 0.60 ± 0.02, which is slightly larger than the value 
for the homopolymer case, meaning that heterogeneity 
is a significant perturbation in this model. For a poly- 
mer with unbalanced numbers of positive and negative 
charges ^ N^) we find that Tg decreases, and even- 
tually vanishes, with increasing the excess charges frac- 
tion X = \ N+ - N_ \ /N (see Figure |l|). Finally, we 
explore the ground state and the energy landscape of a 
neutral (i.e. X — 0) polymer, in an attempt to find evi- 
dence for the existence of a glass-like freezing transition 
for such a polymer. We do not find much evidence for 
this transition, a result which is in accordance with the- 
oretical predictions for our model. 




FIG. 1. Phase diagram of a random heteropolymer in the plane 
of temperature (T) and excess charge (X). Vertical bars indicate 
estimated uncertainties in Tq. The lighter shaded area indicates 
X values where the results are rather ambiguous, due to the poor 
quality of the MC results. 
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Following Kantor and Kardar |l4|, we study poly- 
mers with quenched heterogeneity and a fixed number of 
monomers A'^. The calculation of thermodynamic quan- 
tities is performed as follows: For chain lengths of up 
to 15 steps (16 monomers), we completely enumerate all 
possible spatial conformations, and average over all pos- 
sible charge sequences (quenches) for a given X. The 
number of the spatial conformations grows exponentially 
with N |15|. However, we take advantage of lattice sym- 
metries to reduce the number of independent configura- 
tions. Thus, for a 16-monomer neutral chain we have 
to enumerate 802,075 different configurations (unrelated 
by symmetry), and average further over 12,870 quenches, 
giving ^ 10^" possibilities. For chain lengths of up to 23 
steps, we still enumerate all conformations, but average 
over a limited number (20-100) of quenches. For chains 
of up to = 100, we use MC simulations, applying the 
"pivot" algorithm p6[. 



The presence of a collapse can be clearly seen in Fig. 
H, depicting two conformations of a 50-monomer neutral 
polymer, sampled by the MC procedure at high and low 
temperature. We expect the polymer to be SAW-like 
(with I'sAW = 0.75) at high temperature, and become 
compact (incompact = 1/d = 0.5) when we lower the tem- 
perature. Thus, R^g/N ^ const in the collapsed state and 
R'^/N ~ in the swollen state. Temperature depen- 
dence of R^g/N for various A^s, depicted in Fig. ^, indeed 
exhibits such behavior. The presence of a phase transi- 
tion is also indicated by the peak in the heat capacity, as 
shown in Fig. ^ The peak grows and slightly shifts to- 
ward higher temperatures with increasing N . Out results 
bear qualitative resemblance with the known behavior of 
homopolymers (see, e.g. M). 




FIG. 2. Expanded and collapsed configurations of a neutral 
50-monomer polymer at T = 2 (top) and T = 0.6 (bottom), ob- 
tained by MC simulation. Oppositely charged monomers are de- 
noted by dark and light filled circles. 




FIG. 3. Squared radius of gyration, divided by chain length, 
for neutral polymers. Distances are measured in lattice constants, 
temperature is normalized by monomer interaction energy. Curves 
are for all quenches (unless stated otherwise) and the following 
number of monomers (from bottom right) : 6 , 8 , 10 , 12 , 14 , 16 
, 18 (100 quenches) , 20 (100 quenches enumeration, 50 quenches 
MC) , 22 (100 quenches) , 24 (20 quenches) , 26 (50 quenches), 
36 (25 quenches) , 50 (10 quenches) , 80 (10 quenches) , 100 (10 
quenches). Solid lines represent results of enumeration, connected 
circles represent results of Monte Carlo simulation. 




FIG. 4. Heat capacity divided by chain length, for neutral poly- 
mers. For details of chains see Fig. ^. 

We cannot use the simple method used to find Te in the 
3D case , which is to observe the intersection of the 
graphs of R^/N vs. T for different A^s (using the facts 
that ve = t'idcai = 0.5 and i^sAW > ve > incompact). This 
is because in the 2D case, z^compact — t'ldoai — 0.5. More- 
over, unlike in 3D, it is believed that for 2D, vg > t'idoai 
1^,0, and the 2D value of i^e is not known a priori. The 
method we use to estimate Tg and i/g is based on the 
observation, made by de Gennes, that the 0-point is a 
tricritical point P,p^. Based on the general theory for 
tricritical phase transitions, it has been shown theoret- 
ically |Q,|l9| and verified numerically for homopolymers 
that Rg for different temperatures and chain lengths 
can be described using a single scaling function: 

Rl/N'-' ^ f±{N'f'T) , (2) 

for A^ > 1 and t < 1, where t = |T - Te\/Tg and cj) is 
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a crossover exponent. The scaling function f±{x) should 
have the following limits: 



{const for x ^ 

^2{v^^^~ve)l4> for a; oo , T > Tq (3) 
2,2(>.co„,pact-!/£,)/<;4 for a; ^ oo , T < Te . 



It is easily seen that the asymptotes of f± (x) at high and 
low temperatures reconstitute the behavior of a SAW and 
of a compact chain, respectively. 

Fig. H depicts the scaling function Rg/N'^'^" vs. the 
scaling variable N'^t, using data of all chain lengths (7- 
99) for temperatures satisfying r < 0.25, and choosing 
the parameters Tg , vg and so that the points fall 
on two converging lines as required. As seen in the fig- 
ure, a very good data collapse is achieved for Tg = 0.83 
and vg — 0.60, and the resulting lines also approach the 
slopes of the theoretical asymptotes. It should be fur- 
ther pointed out that the scaling behavior is found to 
be very sensitive to the value of vg, slightly less to the 
value of Tg, and quite insensitive to the choice of (p. Af- 
ter examining numerous figures of the like of Fig. ^, using 
various parameter values, we evaluate the transition pa- 
rameters to be Tg = 0.83 ± 0.02 and vg ^ 0.60 ± 0.02, 
where the errors were estimated according to the param- 
eter values where data did not collapse anymore (accord- 
ing to our subjective judgment). We verify these val- 
ues by plotting Rg/N'^^" vs. T for various chain lengths 
(Fig. . The curves should intersect at T — Tg (because 
'^SAW > vg > ^compact)- This is iudccd what happens, 
within the estimated errors of the parameters and of the 
simulated data. 
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FIG. 5. Logarithmic plot of Rl/N'^^'e vs. N'^'r, close to 
the tricritical temperature. Data from chain lengths 7-99 
is used, with values of t up to 0.25. Parameters used: 
Vg = 0.60, Tg = 0.83, 4) = 0.636. Dots denote enumeration data, 
circles denote MC data. Dashed lines show theoretical asymptotes 
(Eq. p|), with amplitude fitted to data . 
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FIG. 6. B?g/N^'^o vs. temperature, for neutral polymers of 
length 7-99 (For details see Fig. ^). The graphs intersect at the 
9-point. 

Our value of vg seems to be larger than the 2D ho- 
mopolymer value, for which most estimates lie in range 
0.51-0.58 1^^. In order to support this claim, we have 
used our own enumeration and MC procedures (prop- 
erly altered) to simulate homopolymers, and found that 
jyhomo _ Q _|_ Q Therefore, our conclusion is that 
j^hctcro ^ ir^°™° , Simple dimensional arguments show 
that randomness is marginally relevant for ideal polymers 
in 2D. One might expect that swelling of the chain at the 
(?-point will render randomness irrelevant. Beyond sim- 
ple arguments ||lo[] no extensive analysis of the problem 
has been performed. 

What happens to the collapse transition if the charges 
on the chain are not exactly balanced? We have repeated 
the procedures described above for increasing values of 
excess charge X. A decrease in Tg is observed with in- 
creasing X, up to a value X « 0.6 where the collapse 
vanishes and SAW behavior prevails at all temperatures. 
The phase diagram in {X, T) plane (Fig. |l|) is similar to 
that obtained for 3D [H. 

Finally, we have attempted to find evidence of a 
freezing transition for a neutral polymer, analogous to 
the glass-like transition occurring in the random energy 
model, which was observed for other heteropolymer mod- 
els An indicator of the transition ||ll[] is the pa- 
rameter x(T) = 1 — X^fcPfc' where pk is the normalized 
Boltzmann weight of a given conformation k. x{T) char- 
acterizes the number of conformations which are thermo- 
dynamically relevant at a given temperature, and in the 
case of a freezing transition it should decrease - when 
averaged over all quenches - from a; « 1 at the freez- 
ing temperature to x w at T = 0. Investigation 
of our model, however, shows that x{T) does not reach 
zero value even at T = 0, due to the degeneracy of the 
ground state. In addition, the decrease in x occurs in 
the same temperature range where the polymer collapses, 
meaning that the decrease in number of relevant confor- 
mations comes mainly from the folding - as is the case 
of homopolymers. Another attempt to validate the anal- 
ogy with the REM was done by examining the similar- 
ity between the degenerate ground state conformations 
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of a randomly chosen polymer. The similarity between 
each two configurations was characterized by the number 
of monomer pairs which are nearest-neighbors in both of 
the configurations It was found that, on the average, 
there is indeed a significant dissimilarity between ground 
state conformations. However, this dissimilarity is not 
pronounced enough to declare that the ground states are 
definitely "structurally different", as is the case in the 
REM. 

We note, that the seeming absence of a freezing transi- 
tion in our model may come from specific features of this 
model: The lack of a strong attractive homopolymeric 
term in the Hamiltonian, which would lead to the forma- 
tion of a "molten globule" , and a lack of sufficient hetero- 
geneity in the monomer-monomer interaction. These two 
may be required in order to achieve an analogy between 
a heteropolymer and the REM, and are usually included 
in model studies, both analytical and numerical, which 
exhibit a freezing transition |po|| . In addition, it appears 
|f7|,pT| that formation of a unique structure in heteropoly- 
mers is very sensitive to space dimensionality, with d = 2 
being a marginal and nonuniversal case that strongly de- 
pends upon the type of lattice, type of interaction and so 
on. 

In conclusion, we have investigated a 2D lattice model 
of polymers, with a quenched random short-range inter- 
action. We have seen that a neutral polymer under- 
goes a tricritical ^-transition with a critical exponent 
lyg = 0.60 ±0.02, a value that seems to be higher than the 
homopolymer value. Since the difference is rather small, 
further MC studies with larger N, and renormalization- 
group studies, are needed to verify the difference between 
the exponents. For non-neutral polymers, we have ob- 
served a decrease in the ^-temperature with increasing 
excess charge, until the collapse disappears and SAW 
behavior prevails at all temperatures. Finally, we did 
not find evidence for the existence of a glass-like freezing 
transition for a neutral polymer, a result which seems 
to be in accordance with theoretical predictions for this 
model. 
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